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Невзирая на успехи в определений уклонения отвеса по спутниковьїм и наземньтм данньтїюм 
гравиметрии, не решеньт вопросьт разрешимости граничньтх задач для определения превьшшений лока- 
льньх участков геоида по данньгкм силь тяжести. 

Для приближения участков геоида предложена модификация метода Молоденского, которая 
позволяет определять пространственное распределение превьшений геоида, аппроксимируя их точе- 
чньми источниками по данньтм потенциала притяжения, заданньтм в ограниченной области. Опреде- 
леньт соответствующие аналитические конструкции для силь тяжести, потенциала притяжения, ано- 
малий вьгсот геоида. 

Поставлена граничная задача последовательного уточнения рельефа геоида по данньтм потен- 
циала притяжения в ограниченной области, Для еб решения средствами функционального анализа по- 
лучено интегральное уравнение сильт тяжести. Исходя из него, при ряде допущений относительно 
ориентации векторов реальной и нормальной силь тяжести полученьт аналитические аппроксимации 
для описания аномалий вьгсот геоида и уточнения фигурьт Земли - в виде линейного интегрального 
уравнения 1-го рода. 

Исходя из приближений потенциала притяжения, определеньт вьісотьт геоида в виде линейной 
граничной задачи с граничньми условиями в виде нелинейного функционала над банаховьтм прост- 
ранством области наблюдений. Она сведена к линейному интегральному уравнению 1-го рода для по- 
тенциала притяжения - уравнению фигурьт Земли. Предложен итерационньй алгоритм для вьгчисле- 
ния приближений потенциала притяжения. Доказана сходимость последовательньїх итерационньтх 
приближений фигурьі Земли. Зта сходимость считаєтся доказательством существования решения в 
численном виде. Устойчивость вьтчислений обеспечиваєт применение способа Сизикова для числен- 
ного решения интегральньх уравнений 1-го рода. 
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Введениєе. Невзирая на успехи в определениий уклонения отвеса (отклонения поверх- 


ности геоида от сфероида) как по спутниковьм |1|, так и наземньгкм данньтм |2| гравиметриий, 
остаются нерешенньми многиє вопросью теоретического обоснования соответствующих 
методов, и в частности, о корректной постановке и разрешимости граничньтх задач для опре- 
деления превьшений геоида над земньїтм зллипсоидом в различньтх системах отечета. Исходя 
из того, что поверхность геоида под сушей не поддабтся математической формализации, и 
лишь приближенно описьваєтся квазигеоидом Молоденского, и вследствиє динамического 
изменения поверхности геоида после крупньх землетрясений и значительньїх тектонических 
движений, не теряет актуальности разработка новьїх способов определения локальньх участ- 
ков фигурь Земли по данньтм силь тяжести (31. 

Альтиметрические методьт уточнения структурь! геоида (4| хорошо описьгвают участ- 
ки морской поверхности, но имеют ограниченную точность над сушей вследствиє снижения 
точности определения сферических гармоник на вьтсоких степенях разложения в ряд по сфе- 
рических функциях. Кроме того, в модели |1|, полученной зтими методами нет надлежащего 
учета сферически несимметричной части глобального поля сильт тяжести вдоль оси враще- 
ния планеть (ее надлежит учитьвать точечньми источниками). 

Для приближения участков геоида, по крайней мере, его вьсокочастотной области, 
можно применить некоторую модификацию метода Молоденского, которая позволяет опре- 
делять пространственное распределение превьшений геоида, аппроксимируя их точечньшми 
источниками по данньтм потенциала притяжения, заданньм в ограниченной области. Иссле- 
дованиє соответствующих аналитических конструкций предлагаєтся вашему вниманиюї. 
Потенциал притяжения У(х). Пусть на гладкой поверхности СС , разделяющей ограничен- 


ную С" инеограниченную С" части области С. , заданьтг значения сильт тяжести: 


ОМУ(х) 

«(х) з , п(х) 1.08), з Їх: УУ(х)- С). (1) 
дп(х) 

Пусть в области С определен центр тяжести в начале неподвижной системьі декар- 

товьіх прямоугольньх координат Ох, х,х; (ее ограничения определим позже). Если область 


С" вращаєтся вокруг оси Ох, с постоянной угловой скоростью «Фо Ф, то значения «(х) 





силь тяжести на СС равньк: 
в (х) - 2(х)- а (х, со5(п, х, )-к х, соб(п, х,)) хе С. (2) 
Определеннье по (2) значения «(х) будут некоторьтм приближением згау У(х) , Где 
У(х) - М/(х)- О(х), Ух) -- потенциал сильг тяжести, О(х) - центробежньй потенциал, ибо 
нормали п(х) 1 дС,, и п(х)1.дС.,, где ОС, - Їх: У(х)- с), не совпадают. Однако мь! можем 
пренебречь зтим, считая что (п(х), п (х)) -0. 


Ставя задачу - определить по граничньтм данньм (2) потенциал притяжения в виде 
1 
уУ(х)- Ї Ме) 4с,,хе С", (3) 
дл 2. г(х, 2) 
в итоге получаєм для У х є С" уравнениєе сильт тяжести: 


здо а Й по іс пауаа (4) 


На его оснований после ряда функциональньх преобразований получаєм нелинейное 
интегральное уравнение для определения плотности це), с єдС в виде: 


3 ш(х)я а Ї Ї не п іст аа, -(х) хедС. (5) 


4 дос б 
Определив в вьтражений (5) функцию ці), Сє ДС, можем вьїчислить потенциал 


У(х), х є 5(Е2, К) по формуле (3). 
Все вьтчисления можно проводить при условиий, что значения «(х) определеннь. Но, 





2 


мьт не знаєм величин направляющих косинусов сок(п(х), з і - 1,3 внутренней для области 


С" нормали п(х) и, следовательно, не можем определить точно и сами значения «(х) ; 


сох(пу (х), х,) - со5(о, хх. ) ; (6) 
и вьшолним необходимьге вьтшисления. В результате получим приближенньте значения поте- 
нциала притяжения УФХ) ,хеб (2,8). 


Положим на первом зтапе 


2. Аномалий вьгсот геоида Є й В разньгїх точках пространства расстояние между 
зквипотенциальньми поверхностями различно, в соответствии с различньжми значениями си- 
ль тяжести. Для придания физического смьгсла превьшшению одной точки поверхности СС 
над другой еє точкой, например, В над А, или С над А (рис. 1), определим вьсоть в виде 





Рис. 1. Расхождениє векторов реальной и нормальной силь тяжести при определений аномалий вьБісот. 


Но ,2)є РУТ «4, (7) 


где интегрирование "нивелирньх превьшшений" АЙ ведется в поле сильт тяжести вдоль сило- 
вой линий ССС", з, (ф, А) - среднее значение сильт тяжести на отрезке ССС". Так как ве- 


личиньт 8, (фФ, А) мьт не знаєм, вместо вьгражения (7) приходится вьічислять следующеє 
1 

Н, (р.4)з з-ай, (8) 

Й ОЇ 


где уп | ,2) - среднее значение нормальной силь тяжести, т.е., значение у, в точке ( 4) 


при Н -0.5Н ,( ,1). Разность Н, (ф,4)-Н б ,А) будет р аномалию вьгсоть 





Круа)- ню, д)- Н, (руд) - Тк(9.2) - (р.а У «4 (9) 
8 (Ф.)-у, (ФЛ 


Аномалий вьгсот «м ) определяют по методам зи и и Молоденского из 
линейной граничной задачи для уравнения Лапласа с граничньми условиями для искомой 
функции ТМ ), задаваємьтми в виде 


ау) ВОМєт - (м) М 


где М - неизвестная точка земной поверхности с вьшотой Н,, М - точка земной поверхнос- 


З и 1 ду(м) - 
ти с вЬІСОТОЙ Но ам) а (м) ду(м)! в(м)є 1. 





Но указанное соотношение нецелесообразно для практического применения. Его мож- 
но упростить, сделав ряд допущений, в частности, о представлений потенциала силь тяжести 
в виде У(М ) з (М у ТМ ), и со5(п, у) з 1, где (/ (М ) -- нормальньй потенциал, т(М ) - во- 


змущающий потенциал, У(М ) - внутренняя нормаль к уровню потенциала (/ (М ) «С. Тогда 


будем иметь соотношениє 
1 ду(м) от(м) 
- з з2(М)-у(М). 10 
Если бьг поверхность Земли с ортометрическими вьсотами Н, находилась вне гипсо- 
метрической поверхности Х, определяємой нормальньми вьсотами Н,, то решение 
«м )- ТМ у УМ ), М єдС задачи бьло бьг гарантированно точноє. Но для зтого, как сле- 
дует из вькражения (9), нужно «м ) » 0. А зто условие не связано с граничнькм условием (10) 
и для его обеспечения придется менять (10), т.е., решать другую граничную задачу. 
Определим решение исходной задачи (2-5). 
3. Определение рельефа рі(ф, А). Нам известно приближенное значение потенциала 
УФ )- Т б. У; ас. . Введем замену вида 
дл УСоє) 
х, з Кві ф, с05 Л,, х, з К8іпФ,8іп Лу, хз з Ксо5ф,, 
2, - рзіпфсо84, б, - рзіпф5іпЛ, 4, - рео8ф. 
Вьчислим функционал: 
ВОСКУ С ОЮА 333; 
І (р, ф,А)- з з 
д(р.1)Ї 1д(9.4)1 1 дію) | 


Найдя соответствующие частнье производньєе, получим 
2 2 2 2 
Я др) . др др б др 
І"(р,ф,А)- р о віп'ф-| 3 | віп (22) -р зерно іо 6 ей М С 
(р,ф,2)- р ї екзо З Де НеНН оч за 


2 
и, следовательно, ас. «Де м гч(Р) адал, дб рал,О«АкОлу. 























дФ 


Таким образом, задача (2-5) получила граничнье условия в виде нелинейного функци- 
онала над пространством С (аас), виь и дифференцируемого при р «К: 











др д 
УВУ род мае бо Фу рід, (1) 
| 2 2 
| д д 
др д моли Р угон 28) 
где Дод, 2) А вна , рері(т,4), 
од 4К'яр'-2Кроозу 


СО5 у з С05 ФСО5 Ф, -5іп Ф5іп Ф, со5(2 - Ло). 
После ряда аналитических преобразований получим разность потенциала притяжения 
1 д 0? д 
АУ(у.р)- | а оо Ума 
АлояА00  ддби  дАбу 


Отсюда легко получить линейное интегральное уравнение 1-го рода для искомой функции 


Гкбоє, р, (р, АМмема з т(х), х є (2, В), (12) 
дГ(р,) "То, ) Ф"Т( 5 | 


др | дфіи || даду 








1 
где тер,я -Ри» к(х,2,)е І 


бачити), яке а | Коба 


ди ди Орди др ди 1 би, 1 
др дф др дА др  дф др дл др 
дФ ОА 
Вьчислив У" Р 9 ) і -- 0,1,2,..., определим по формуле (12) значение рід (ф, А) и затем в 


пчі 
точках зтой поверхности 9с - направляющие косинусьт сов(л(х), хі), ) -1,2,3. Далее 
зро следующеє приближение граничного условия (2) для отьюскания потенциала 
ув 9) і - 04,2,.. .ит.д. 


Вводя обозначения: 


Ран -рт,, МУ, (ру, и, - др, 2 У чі др, 7 а з (Ку фо, дь)є 5, В -(р.ф.4)єбс, 








ї 2 2). :12 2; 
збору орелуо«а«зя), нем обоє азноо 1 мор зи )зіп Фу (3) 
Дл «К'юр'-2Кроозу 








Со5 у - с05 ФСО5Ф, -5іп Ф5іПФ, со8(Л - Ло; 
получаєм следующее уравнение фигурь Земли 


Г ку (с, Вт (РІВ - У(а)- а Гео, ув, 14) 


ее 





ов 1 |2к(с УУ,) О'Е(с, У) д'К(а, М) 

з 4л др дфди ОАбм 

4. Сходимость метода. Доказана сходимость последовательньх приближений (14) 
фигурь Земли, представив его в виде последовательного итерационного уточнения: 


Тк а, Ву (РАВ - Ї ка, 8) ть у, Па )- (о. (15) 


Обозначив зкстремумь ядра К, (а, В) как тіп 





Гк (а,виВ 





» Ко и тах Гк/а вай «К,, 
дос 


будем иметь для оценки искомой разности цепочку неравенств: 


рід) р), «(8 рр зо) інею 


Кк Кк 


откуда и получаєм в силу условий (12), (13) для вьшуклой поверхности ОС доказательство 
сходимости уравнения (14). 

Вьводь. Исходя из линейного интегрального уравнения силь тяжести (4), представ- 
ленного в виде нелинейного уравнения для восстановления плотности, полученьт аналитиче- 
ские аппроксимации для описания аномалий вьсот геоида (9)-(10) и для уточнения фигурьі 
Земли - в виде линейного интегрального уравнения 1-го рода (14). Для последнего указан 
итерационньй способ нахождения последовательньх приближений (15) и доказана его схо- 
димость (18). Последнее обстоятельство при исходньх предположениях о гладкости границь 
области СС. и ограниченности ядра К, (а, В ) следует считать и доказательством существова- 





ния решения в численном виде, исходя из предложения Маслова |5). 

Ради обеспечения численной устойчивости алгоритма (14) следует избрать для его ди- 
скретизации и вьшчислений подход, изложенньї(й в |6| для численного решения интегральньх 
уравнений 1-го рода, в частности, заменив интегральное ядро К, (а, В ) интегральной суммой 


по методу трапеций при разбиении интервала интегрирования равномерной сетью точек. 
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